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Spectroscopic Search for Faint Tertiaries 
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ABSTRACT 

It is unclear how very close binary stars form, given that during the pre-main-sequence phase 
the component stars would have been inside each other. One hypothesis is that they formed 
further apart, but were brought in closer after formation by gravitational interaction with a third 
member of the system. If so, all close binaries should be members of triple (or higher-order) 
systems. As a test of this prediction, we present a search for the signature of third components in 
archival spectra of close binaries. In our sample of 75 objects, 23 show evidence for the presence 
of a third component, down to a detection limit of tertiary flux contributions of about 0.8% 
at 5200 A (considering only contact and semi-detached binaries, we find 20 out of 66). In a 
homogeneous subset of 58 contact binaries, we are fairly confident that the 15 tertiaries we have 
detected are all tertiaries present with mass ratios 0.28 < M3/M12 < 0.75 and implied outer 
periods P < 10^ d. We find that if the frequency of tertiaries was the same as that of binary 
companions to solar-type stars, we would have expected to detect about 12 tertiaries. In contrast, 
if all contact binaries were in triple systems, one would expect about 20. Thus, our results are 
not conclusive, but sufficiently suggestive to warrant further studies. 

Subject headings: stellar dynamics — methods: data analysis — binaries: close — stars: formation 



1. INTRODUCTION 

Most stars are in binaries, yet our understand- 
ing of how these form is far from complete (for a 
recent review, see Tohline 2002). Particularly puz- 
zling is the existence of very close binaries, with 
orbital separations of just a few stellar radii. That 
these cannot form independently is easily seen: 
during the pre-main-sequence phase, these stars 
would have been inside each other. Yet they exist. 

One hypothesis is that they were originally sin- 
gle entities which spun up to break-up velocity 
during contraction and split in two. It is un- 
clear, however, whether this world work: descend- 
ing the Hayashi track, stars are centrally concen- 
trated and fission into roughly equal parts appears 
implausible. 

If close binaries cannot form directly, could 
the stars perhaps form as wider binaries, and be 



brought closer together later? One such possi- 
bility is that the binary is part of a hierarchical 
triple, and shrink due to interaction with the third 
component (Kiseleva 1998; Eggleton & Kiseleva- 
Eggleton 2001). This can work as follows: the 
tertiary induces Kozai cycles (Kozai 1962) in the 
inner binary, in which angular momentum is trans- 
ferred between the inner and outer system, leading 
to cycles in eccentricity and relative inclination. 
For point masses this process is cyclical, but for 
stars, if the eccentricity becomes sufficiently high, 
tidal effects take over at periastron, and the or- 
bit will circularize with a final separation of about 
twice the periastron distance, i.e., much smaller 
than the initial one. 

What helps the above is that Kozai process is 
weak: it only works in the absence of anything 
else. Thus, no cycles occur while there is tidal 
interaction between the stars and/or their disks. 
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i.e., as long as the stars arc young and big, nor 
once they have been brought in close together. 
The only requirements are that a sufficient num- 
ber of binaries be members of triple systems, and 
that many have sufficiently small initial inner sep- 
arations and sufficiently large relative inclinations 
between the inner and outer binary planes. None 
of these constraints appear problematic: 15-25% 
of all stellar systems have three or more com- 
ponents (Tokovinin 2004), and some well-known 
triples have high relative inclination (e.g., Algol; 
Lestrade et al. 1993). 

Indeed the mechanism has been invoked to ex- 
plain the properties of a number of individual sys- 
tems, such as the triple TY CrA (Beust et al. 
1997) and the quadruple 41 Dra (Tokovinin et al. 
2003). Furthermore, it was found by Tokovinin & 
Smekhov (2002) that many visual multiples had 
close spectroscopic subsystems, which might have 
formed by the above process. 

The hypothesis has not, however, been taken 
to its logical conclusion: do all close binaries form 
this way? The beauty of this perhaps far-fetched 
suggestion is that it makes a very clear prediction: 
all close binaries should be in hierarchical triples 
(or higher order systems). In this paper, we inves- 
tigate this possibility with a particular subset of 
close binaries, the W UMa contact binaries. 

W UMa contact binaries - in which the two 
companions share an outer envelope - are the clos- 
est known binaries. While these have had yet a 
further phase of orbital shrinkage, likely related to 
magnetic braking and/or gravitational radiation, 
this phase could only happen if they were very 
close binaries to start with (Vilhu 1982; Stepien 
1995, and references therein). Intriguingly, many 
contact binaries appear to be accompanied by ter- 
tiaries: for instance, Rucinski & Kaluzny (1982) 
noted the frequent presence of visual companions, 
while Hendry & Mochnacki (1998) found the spec- 
tral signature of a tertiary in a number of contact 
binaries. Furthermore, in radial-velocity studies 
of contact and other close binaries aimed at mea- 
suring their orbital parameters, one of us (S.M.R.) 
has found that about one in four binaries showed 
the signature of a tertiary component in its spec- 
trum (see Rucinski 2002 and other papers in the 
same series). 

The above led PribuUa & Rucinski (2005, here- 
after Paper I) to collect the available evidence for 



multiplicity for contact binaries with V < 10. For 
the better-observed Northern-sky subsample, they 
inferred a multiple frequency of 59 ± 8%. Since no 
method can detect all miiltiplcs, this is a lower 
limit to the true fraction, but it is difficult if not 
impossible to extrapolate given the complex selec- 
tion effects for various techniques and companions 
types. 

Here we present a detailed analysis for one 
particular technique of searching for the spectral 
signature of tertiaries. We re-analyze the data 
sets used for the radial- velocity studies referred to 
above using a new technique, outlined in Sect. 2, 
optimized for the detection of tertiaries. We dis- 
cuss possible pitfalls and systematic effects, and 
find that these limit us somewhat, but that we 
can still detect tertiaries down to flux ratios of 
about 1%, an improvement by a factor of three 
compared to the earlier results. In Sect. 3, we infer 
properties of tertiaries and check consistency with 
previous work. We discuss limits and biases in 
our sample in Sect. 4, and use the companion dis- 
tribution of solar-type stars measured by Duquen- 
noy & Mayor (1991) to test the hypothesis that all 
close binaries are in triple systems. In Sect. 4, we 
also discuss what the alternative, null hypothesis 
should be, i.e., what one would expect if multiplic- 
ity had little or no influence on the formation of 
close binaries; conservatively, we assume that for 
this case the tertiary frequency be similar to the 
companion frequency of regular stars. We sum- 
marize our results and discuss future prospects in 
Sect. 5. 

2. DATA SET AND ANALYSIS TECH- 
NIQUE 

The data set we have available is that used in 
papers I-IX of a series on "Radial Velocity Studies 
of Close Binaries" (for an overview, see Rucinski 
2002). Briefly, it consists of almost 4000 spec- 
tra of 75 close binaries, all taken at the David 
Dunlap Observatory. We list the objects in Ta- 
ble 1. For our purposes, the important character- 
istics of the spectra are that they were taken with 
the 1800 line mm^^ grating centred at 5180 A (cov- 
ering about 200 A around the Mg I triplet), and 
through a 1'.'5 or 1'.'8 slit (matched to the typical 
seeing of 1'.'7). The resulting slit images project to 
0.64 and 0.80 A, respectively, and, including the ef- 
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feet of seeing, lead to an effective resolutions of 35 
to 50knis~^. 

We search for the spectroscopic signature of a 
tertiary by fitting the spectrum of the contact bi- 
nary, and checking whether adding a spectrum 
of a fainter tertiary improves the fit significantly. 
For the binary spectrum, we can make use of the 
convenient fact that the contact binary has a sin- 
gle spectral type (due to energy transfer from the 
more massive to the less massive component, by 
a mechanism not entirely understood; for a re- 
cent review, Webbink 2003), and that its lines are 
strongly broadened by the rapid rotation and or- 
bital motion. In contrast, the lines of the third star 
should be narrower since there is nothing to have 
prevented it, like all low- mass stars, from slowing 
down. (Note that there are exceptions: new DDO 
observations have revealed a broad-lined A-type 
companion for V752 Mon, which entirely masks 
the radial velocity signatures of the binary so that 
only the variability is detectable. Since we remove 
more massive, early- type stars, this will not bias 
our sample.) Furthermore, any orbital motion of 
the tertiary should be small compared to that of 
the contact binary, which implies that we can ana- 
lyze spectra averaged over the orbital phase of the 
contact binary. This not only increases the signal- 
to-noise ratio, but should also smooth out possible 
relatively sharp- lined features from the contact bi- 
nary, such as might be caused by star spots. 

Below we describe how we implemented the 
technique and how we model the line broadening 
in the contact binary. We then discuss the crite- 
ria we use to determine whether a detection of a 
tertiary is significant and real, and determine our 
sensitivity limits. We conclude with a discussion 
of the limitations of our method. 



der to find the one that matches best. The model 
broadening function is composed of three Gaus- 
sians with equal width. We choose such a model 
to ensure that our broadening function is wide 
enough that the sharp-lined signal from a putative 
tertiary is not removed; since this was a particu- 
larly difficult part of our analysis, we describe it 
in more detail in Sect. 2.2 below. 

To measure the quality of the fit, we convolve 
the model spectrum with a truncated Gaussian (to 
simulate the effects of seeing and transfer through 
the slit), regrid on the observed pixel array, mul- 
tiply with a polynomial function to simulate dif- 
ferences in the normalization, and finally calcu- 
late the variance between the observations and 
model. We minimize the variance as a function 
of the various parameters using the Downhill Sim- 
plex method as described by Press et al. (1992, 
§10.4). 

Once the best fit to the contact binary spec- 
trum has been determined, we add a third star to 
the spectrum, optimize the relative flux and ve- 
locity, and determine the resulting improvement 
in the fit. We again repeat this procedure for a 
wide range of different spectral types (from early 
F down to early M). 

With our procedure, the final set of parameters 
determined is: (i) best-fit template spectrum for 
the contact binary, or, more interestingly, its tem- 
perature; (ii) the systemic velocity; (iii) the width, 
separation, and two relative intensities for the 
three-Gaussian broadening profile (see Sect. 2.2); 
(iv) the best-fit template for the third star; (v) its 
fractional intensity; and (vi) its velocity relative 
to the contact binary. In addition, there arc up 
to 10 parameters without physical meaning, viz., 
those that describe the polynomial accounting for 
difference in continuum normalization. 



2.1. Implementation 

To search for tertiaries, we fit the data us- ^.2. The Broadening Profile 

ing the procedure outlined in Fig. 1. Starting We found that the most diflacult part of our 
with the average of all normalized, barycentered analysis was to accurately reproduce the line- 
spectra of a given source, we try to reproduce it profile shape of the contact binary In princi- 
with a template spectrum of a sharp-lined, slowly pie, with good phase coverage, one might expect 
rotating star convolved with a model broaden- that a single Gaussian would sufiice. In practice, 
ing function optimized to match the binary's line however, the binaries were observed preferentially 
profile. For our templates, we use the database near the quadratures. The result, for systems with 
of high-resolution, normalized stellar spectra ob- ^lass ratio near unity is a double-humped line pro- 
tained with ELODIE (Prugniel & Soubiran 2001), gig in the average spectrum, while for systems with 
repeating our fit procedure for all spectra in or- 
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Fig. 1. — Search for a tertiary component in the spectra of the W UMa system CK Boo. Panel A, from top 
to bottom: /obs, the observed, average spectrum; /12 a template of similar spectral type used to represent the 
contact binary; /obs, /12 (81 BFh: the observed spectrum (thin line) overdrawn with the template convolved 
with a best-fit broadening function (thick line); /obs, /12 ® BFt + fs: the observed spectrum with overdrawn 
the best-fit model composed of the template convolved with a re-fit broadening function, plus a tertiary 
spectrum; f^: the best-fit tertiary contribution. Panel B: Comparison between the residuals from fitting 
the observed spectrum without including a third star (offset by a constant value), and the best-fit tertiary 
spectrum. Panel C: The broadening function used to represent the line profiles of the contact binary. Crosses 
indicate the empirical broadening function found from least-squares decomposition (see Rucinski 2002), and 
the dashed line the fit to those points with our three-Gaussian model shape. This fit is used as an initial 
guess for line profile; the dotted line represents the final shape, after convergence of our procedure. Panel 
D: variance of the fit residuals as a function of the velocity of the tertiary. For CK Boo, this shows only 
one minimum, which is close to the systemic velocity of the contact binary, as expected for a real tertiary. 
Panels E and F: variance of the residuals as a function of temperature of the template used for the contact 
binary and the tertiary, respectively. For CK Boo, the4best fit tertiary is substantially cooler than the binary, 
as expected if it is much fainter. Panel G: the best fit relative velocity between the tertiary and contact 
binary, as a function of tertiary temperature. For CK Boo, this is close to zero, as expected for a physically 
associated component. 
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Fig. 2. — Like figure 1, but for the W UMa system EF Boo (fobs), in which we do not detect a third star 
down to a level of 0.8 %. 
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extreme mass ratios the profile shows a prominent 
central hump (from the more massive component), 
as well as two 'side lobes' (from the less massive 
one). 

Based on this structure, we chose to model the 
broadening profile with a set of three Gaussians, 
but constrained to have identical width and sep- 
aration. This leaves four free parameters: the 
width, the separation, and the relative intensities 
of the two outer Gaussians (since the whole profile 
is normalized, the intensity of the central Gaussian 
does not have to be specified). 

Another problem that arose in automating our 
procedure was that the profiles of the different sys- 
tems were so varied that it was impossible to use 
a single set of initial guesses for the broadening 
profile that worked for all systems. To circum- 
vent this, we first determined empirical broaden- 
ing profiles at the instrumental resolution by least- 
squares decomposition (using the technique de- 
scribed in detail by Rucinski (2002); see Fig. 1, 
panel C). Next, we fit these empirical broadening 
profiles using our three-Gaussian model, and used 
the resulting parameters as initial guesses for our 
main procedure. 

2.3. Detections and Detection Limits 

With the fit results in hand, we need to deter- 
mine whether or not a possible improvement re- 
sulting from adding a tertiary is significant. For 
some objects, this is trivial: the residuals from the 
fit with the binary model show a clear signature of 
a different spectral type. Those, however, would 
typically have been found already in the earlier 
studies, since the tertiary component would lead 
to a narrow peak in the empirical broadening func- 
tion (see Rucinski 2002). 

For fainter tertiaries, one could in principle use 

statistical tests to determine whether the improve- 
ment in variance (or, equivalently, x^) is signifi- 
cant. This only works, however, if the variance is 
dominated by measurement noise. In practice this 
is not the case: the quality of the fit is usually 
limited by a mismatch between our model and the 
true contact binary's spectrum. Indeed, even for 
our spectra with the worst signal-to-noise (such as 
UX Eri; Fig. 3, bottom panel), we find that sys- 
tematics dominate. As a result, the quality of the 
fits is not good in a statistical sense, and the use 



of becomes meaningless. 

In Fig. 3, we show the types of more severe sys- 
tematic mismatches that limit our sensitivity to 
tertiaries. The first is a poor match to the binary's 
spectral type, which happens mostly for cooler 
temperatures for which the ELODIE archive con- 
tains relatively few suitable templates. As can be 
seen for the case of AH Aur in Fig. 3 (top panel), 
the mismatch leads to low- frequency residuals and 
thus an increased variance. Since the residuals 
have long wavelength, one could still tease out the 
signature of a tertiary, but since the limits will not 
be as good, wo decided not to include objects like 
AH Aur in our statistical sample. 

The second source for systematic error is more 
problematic: poor matches to the binary star's 
line profile. For cases such as V351 Peg (Fig. 3, 
middle panel), where our three-Gaussian broad- 
ening function does not match the intrinsic profile 
very well, high-frequency residuals are left. In con- 
sequence, there is an obvious danger of a false de- 
tection of a "third star" that matches these resid- 
uals. 

In order to avoid the above pitfalls, we decide 
whether or not a detection is significant using not 
only visual inspection of the residuals, but also 
the following two physical arguments. First, for 
a faint tertiary, the temperature should be sub- 
stantially lower than that of the contact binary, 
and, therefore, as a function of tertiary tempera- 
ture, minimum variance should occur at low val- 
ues. Second, the orbital motion of the tertiary 
should be relatively small, and hence as function 
of relative velocity, minimum variance should oc- 
cur near zero. In the example shown in Fig. 1, 
both criteria are met, and hence we consider the 
detection of the tertiary secure. There are also, 
however, a fair number of sources for which ter- 
tiary flux and the decrease in variance are similar, 
but for which the inferred tertiary temperature is 
higher than that of the contact binary and/or the 
radial velocity is inconsistent. 

With the above procedure, we find that we are 
able to detect tertiaries down to fluxes of 0.8% of 
that of the contact binary. We confirmed this by 
adding third stars at different flux levels to objects 
for which we did not detect tertiaries, and finding 
the level at which we could recover those: we found 
a limit of 0.9% even for cases where the match to 
the line profile was relatively poor, such as UX Eri 
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Fig. 3. — Three examples of problems encountered in fitting the contact binary, which make detecting a 
third star more difficult. Top: AH Aur and its best fit (in bold overlay), with the residuals (enlarged) below. 
A poor match in spectral type to the contact binary leads to large-scale residuals. Middle: V351 Peg and 
the residuals from its best fit. In this case, the dominant residuals come from errors in fitting the contact 
binary's complicated broadening profile, and so show up as high frequency residuals. Bottom: UX Eri and 
the residuals from its best fit. In this case, the spectrum itself is noisier (in addition to the poor spectral 
match) and so the residuals are larger. In general, even with these effects we can detect a third star down 
to a level of about (3 ~ 0.8%; for the worst cases (which we exclude in our analysis) this level is closer to 
/3 ~ 2 - 3%. 
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and V351 Peg. 

2.4. Limitations 

Apart from the problems addressed above in 
obtaining an adequate fit to the binary, our tech- 
nique also has limitations inherent in the assump- 
tions we make about the tertiary. An obvious one 
is that we cannot detect compact objects, since 
these would be too faint and likely not contribute 
any spectral features. Another is that we have 
few templates for late-type stars. This is not an 
issue for the more massive, earlier contact bina- 
ries, for which such late-type tertiaries would be 
undetectable. But for later-type contact binaries, 
we might miss cool tertiaries or, more likely, over- 
estimate their temperatures and fluxes (the latter 
since the strength of the band heads, etc., gener- 
ally increases with decreasing temperature). 

A different limitation arises from our assump- 
tion that the tertiary rotates slowly. For early- 
type tertiaries, this may not be correct. Those, 
however, would be very bright and hence noted 
independently (furthermore, we will exclude them 
from our sample since we cannot be sure the sam- 
ple of contact binaries with such bright tertiaries is 
complete; sec Sect. 4.1). For late-type tertiaries, 
slow rotation is expected unless the star is in a 
close binary and is kept corotating by tidal forces. 
Thus, our procedure will not identify close bina- 
ries as companions (such as the quadruple system 
composed of two contact binaries, BV Dra and 
BW Dra; Rucinski & Kaluzny 1982). In order 
for the projected rotation velocity to be below our 
resolution, i.e., vsini < 50kms~^, one requires 
-frot ^ 1 d (for a 1 Rq star) . But at such short 
orbital periods, orbital velocities would be even 
higher, and those would smear the signal as well 
(at least for systems for which the spectra were ob- 
tained over an extended period of time). Orbital 
velocities for the tertiary were indeed found in HT 
Vir (Lu et al. 2001). In order not to decrease our 
sensitivity, a tertiary that is itself a binary should 
have a radial- velocity amplitude K < 50kms~^, 
which requires Porb ^ 20d (for two IMq stars). 

3. RESULTS 

Out of 75 systems, we have detected tertiaries 
for 23, nine of which had been missed in the orig- 
inal analysis of the spectra. All detections are in- 



dicated in Table 1 and described in more detail in 
the Appendix. Below, we will compare our results 
to those in the literature, and then proceed to in- 
fer tertiary masses and mass ratios. A summary 
of observed and inferred properties of the triple 
systems is given in Table 2. 

3.1. Comparison to Previous Results 

In Table 1, we indicate for all stars in our sam- 
ple whether it was also studied in Paper I, and, if 
so, whether there was independent evidence for it 
being member of a triple (or higher-order) system. 
For many of our detections there is independent 
evidence for multiplicity, and, conversely, for most 
systems for which we detect no tertiary, there is lit- 
tle evidence to the contrary. This likely reflects the 
fact that most methods, whether detecting com- 
panions through gravity or flux, require relatively 
similar minimum masses. 

There are nine systems, however, for which a 
tertiary was found in Paper I but not in our anal- 
ysis. For five cases - EE Get, V2150 Cyg, QW 
Gem, AQ Psc, and AH Aur the discrepancy is 
simply that the separation is too large for any light 
of the tertiary to have entered the slit (for EE 
Get and V2150 Gyg, some light did enter the slit, 
which, knowing that they were visual binaries, we 
could detect; see the Appendix). For three others 
- AB And, V523 Gas, and UX Eri - the minimum 
mass inferred from the arrival-time variations im- 
plies a flux below our threshold (which is the case 
for A'h/McB < 0.28 [Sect. 3.2]; V523 Gas also was 
fit particularly poorly [Sect. 4.1]). For the remain- 
ing system - DN Gam - the identification is based 
on suspected multiplicity from Hipparcos and X- 
ray emission in excess of expectations. Since this 
does not yield a mass estimate, we cannot check 
whether our non-detection makes sense. 

Turning now to the properties of the systems in 
which we detect tertiaries (Table 2), we see that 
the contact-binary temperatures inferred from our 
spectral fits (column Tfit) are generally in fair 
agreement with those inferred from B~V {Tb-v, 
found by interpolation in Table 15.7 of Gox 2000). 
This gives confidence in our method. There are 
three exceptions: KR Gom, 44 Boo and V899 Her, 
in both of which the tertiary contributes a signif- 
icant fraction to the system's light. As a result, 
an incorrect temperature for the contact binary 
is found in the first step of our procedure. Since 
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VIII 




N 


SX Crv 


\" 






S\V Lyu 


1\" 


0.1!)i 




XS\" 22:-!'- 


will 







^13 is the flux ratio (= /3//12) determined from our fit (/3 < 0.008 for non-detections; further properties for the detected systems 
arc listed in Table 2. A brief description can be found in the Appendix for all systems with detections as well as some interesting 
triple systems we missed, marked with f.) 

''Notes *:the spectroscopic signature of the tertiary was not recognized in the original DDO series paper; Y, •, N: the system 
is in the sample of Paper I. parentheses indicate the extended, Vmax > 10 sample: Y-identified as a triple independently from the 
DDO spectra,*— identified as a triple based on the DDO spectra, N— not identified as a triple; NC, BF, TB: the system was not 
included in our statistical analysis because it is not in contact (i.e., the temperatures of the two compoiicuts differ), because our 
procedure yielded a bad fit to its spectrum and hence our sensitivity to tertiarics was poor, or because the tertiary is too bright 
to be sure the sample is complete; TW: the system is included in our statistical analysis, but the detection of the tertiary is not 
counted because the tertiary is at too wide a separation, and we cannot be sure we could detect such tertiaries for all objects in 
our sample. 

"^Variable names: 44 Boo=i Boo; HD 93917=VY Sex; NSV 223=DZ Psc 

References. — I: Lu & Rucinski (1999); II: Rucinski & Lu (1999); III: Rucinski et al. (2000); IV: Lu et al. (2001); V: Rucinski 
et al. (2001); VI: Rucinski et al. (2002); VIII: Rucinski et al. (2003); IX: Pych et al. (2004). 
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this temperature is not checked a posteriori, it is 
expected that the other results will also be inaccu- 
rate (we do not use these systems in our statistical 
analysis; see Sect. 4.1). 

Comparing the flux ratios /3fit = /3//12 from 
our fit to those from the literature (Ait.), there 
are also a number of discrepancies. For resolved 
systems, the literature values should be reliable, 
and hence we need to understand what went wrong 
in our procedure. We see two possible causes for 
errors. First, as above, for systems with bright 
tertiaries (44 Boo and KR Com) the temperature 
assigncid to the contact binary will be incorrect 
and hence the other parameters will be inaccurate. 
Second, for systems with wide separations (V776 
Cas, KP Peg, and, to a lesser extent, 44 Boo), the 
tertiary would have been only partially in the slit, 
and hence the flux will be underestimated. 

Turning now to unresolved triple systems, 
which were all identified in the DDO program 
and thus have values based on the same data, we 
find that our flux ratios are systematically smaller, 
especially for fainter tertiaries. In the DDO pro- 
gram, the tertiaries were recognized by the appear- 
ance of a sharp feature in the broadening function 
at the system's average radial velocity, and the 
flux ratio was determined from the ratio of the 
area under the sharp peak to the remainder of the 
broadening function (Rucinski et al. 2002). Since 
these broadening functions are derived using a 
least-square decomposition based on a single tem- 
plate spectrum, the contribution of the tertiary 
was effectively measured under the assumption 
that it had the same spectral type as the contact 
binary. If its true spectrum has stronger lines - 
as will be the case for faint tertiaries with cooler 
temperatures - this will lead to an overestimate 
of its contribution. Since our procedure uses a 
separate spectral type to derive the contribution 
from the tertiary, our flux ratios should be more 
reliable. 

The above issues allow one to understand the 
discrepancies between literature values and those 
derived here, but make it difficult to estimate re- 
liable uncertainties. From the comparison with 
resolved systems, uncertainties of < 15% are in- 
dicated for flux ratios between 0.05 and 0.5, but 
errors increase towards higher and lower values. 
For the brighter tertiaries, we will use flux ratios 
from resolved observations, which should be good 



to < 10%. For fainter ones, however, one needs 
to keep in mind that our uncertainties increase 
rapidly, reaching ~50% for ratios below ~0.02. 

3.2. Inferred Properties 

Contact binaries follow a period-luminosity- 
color relation, which allows one to derive the ab- 

sohitc magnitude My from the period P (in days) 
and dereddened color {B — V)o (Rucinski 2004, 
and references therein), 

Mv = -AAAlogP + '3m{B -V)o + 0.12. (1) 

With the tertiary flux ratio, this yields the ab- 
solute magnitude of the tertiary. Next, we use 
the fact that contact binaries are on the main se- 
quence, and that, therefore, the generally fainter 
tertiary should be on the main sequence as well. 
Then, with the main-sequence mass-luminosity re- 
lation (we use Cox 2000, Tables 15.7 and 15.8), the 
tertiary mass M3 follows from the absolute V-band 
magnitude. The results of this procedure arc listed 
in Table 2. Here, we did not correct for the typi- 
cally very small reddening [Eb-v = 0.00 . . . 0.03; 
Rucinski, unpublished work). The resulting er- 
rors in the absolute magnitude inferred from Eq. 1 
are on the order of 0.1 mag, substantially below 
the 0.25 mag scatter in the period-luminosity rela- 
tion.^ 

The uncertainty in the derived masses has con- 
tributions from all steps. The magnitudes pre- 
dicted from the period-luminosity-color relation 
are uncertain by ~ 0.25 mag (Rucinski 2004). For 
relatively bright tertiaries, the uncertainty in the 
tertiary flux ratio is smaller (< 15% , or 0.15 mag). 
With a total uncertainty in My of ^ 0.3 mag, this 
leads to an uncertainty in the derived masses of 
< 10%. Evolution along the main sequence will 
likely contribute less, except for the brightest ter- 
tiaries. For fainter tertiaries, the uncertainty in 
the flux ratio leads to an error in My of as much as 
0.7 mag, but since the mass-luminosity relation be- 
comes much steeper for fainter objects, we still ex- 
pect the flnal uncertainty in the mass to be around 
-10%. 

^Ignoring reddening leads to a small systematic underesti- 
mate of the brightness and thus the mass of the tertiary. 
Since the binary's mass will be underestimated as well, 
however, the effect on the mass ratio should be small. 
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Table 2 

Observed and Inferred Properties for Close Binaries with Tertiaries. 



Contact Binary Tertiary 



kJtcU 


P 
(d) 


Ma /Ml 




1 — y 


-'fit 


I'b—v 


My 


Ml 
(Mo) 


Ait. 


Sep. 

("\ 


^fit 


-ffit 


^0 




Ms 


Ms/Ml 


M3/M12 


GZ And 


U. 


.oUo 


U. 


A 

.014 


0, 


.79 


5600 


5200 


A 
4. 


.<5U 


U. 


Oft 


u. 


.Uiilf 


2.13 


U. 


.UiO 


3500 


4000 


y. 


.oy 


n 
U. 


/IO 

.4z 


n 

U. 


A Q 
.4o 


n 
U. 


on 

.zy 


HV Aqr 


U 




U. 


'\ AK 
.140 


0, 


.63 


5800 


5800 


q 
O. 


oo 

.yz 


±, 


.zz 








U. 


noo 
.Uzz 


4000 


4000 


o 
O. 


.UO 


n 
U. 


cn 

.oy 


n 

U. 


A C 


n 
U. 


/IO 

.4z 


V410 Aur 


U. 


.OOD 


U. 


^AA 
.144 





.56 


5400 


6000 


q 
O. 


. / o 


1 

±, 


on 

.zy 


U. 


.ooj 


1.7 


U. 


.zz 


5200 


5600 


A 

4. 


.c5U 


n 
U. 


.y t 


n 

U. 


'TK 
. (0 


n 
U. 


{id 

.00 


44 Boo 





.268 


0. 


.487 


0, 


.94 


5400'' 


4800 


5. 


.50 


0. 


.85 


2. 


.08t 


1.7 


0. 


.23*''^ 


6100'' 


5953 


4. 


.70 


0. 


.99 


1, 


.17 


0. 


.79 


CK Boo 





.355 


0. 


.111 


0, 


.54 


6600 


6100 


3. 


.75 


1. 


.29 


0. 


.007| 


0.12 


0. 


.009 


3900 


3900 


9. 


.14 


0. 


.47 





.37 


0. 


.33 


FI Boo 





.390 


0. 


.372 





.64 


5800 


5700 


3. 


.o7 


1 


.24 








0. 


.012 


3900 


3800 


o 

o. 


.07 


0. 


ro 

.52 





A O 

.4z 


0. 


.31 


AO Cam 


U 


.ooU 


U. 


/lie; 
4io 


0, 


.58 


5800 


5900 


A 

4, 


.Ui 


1 


1 Q 








U. 


UUo 


4200 


4000 




Z-J 


n 
U 


,40 


n 
U 


QO 


U. 


OQ 

.zo 


V776 Cas 


U 




U. 




0, 


.47 


6500 


6500 


Q 
O. 


1 o 


1 
1 


A O 


U. 




5.38 


U. 


UiO 


6100 


5600 


A 


Do 


1 

1, 


,UU 


U 


.0 / 


U. 


.oy 


V445 Ccp 


0, 


.449 


0. 


167 


0, 


.12 


7400 


8400 


2, 


.03 


1 


.95 








0. 


,055 


5600 


6600 


5. 


18 


0, 


,91 


0, 


.46 


0, 


.40 


KR Com 


0, 


.408 


0. 


,091 


0, 


.52 


6100 


6200 


3, 


.42 


1, 


.42 


0, 


.58t 


0.119 


0. 


,23 


6100 


5800 


4. 


,01 


1, 


,19 


0, 


.84 


0, 


.77 


V401 Cyg 


0, 


.582 


0. 


,290 


0, 


.3 


6700 


7300 


2, 


,07 


1, 


.93 


0, 


.03 




0. 


,015 


4700 


4700 


6. 


,63 


0, 


,73 


0, 


.38 


0, 


.29 


V2082 Cyg 


0, 


.714 


0. 


,238 


0, 


.31 


7000 


7200 


1, 


,71 


2, 


.14 








0. 


,020 


5100 


5200 


5. 


,95 


0, 


,79 


0, 


.37 


0, 


.30 


V899 Her 


0, 


.421 


0. 


,566 


0, 


.48 


6300*^ 


6400 


3, 


,24 


1, 


.44 


1, 


.5 




0. 


,73*^ 


6400'' 


6500 


3. 


,59 


1, 


,36 


0, 


.94 


0, 


.60 


ET Leo 


0, 


.347 


0. 


,342 





.61 


5800 


5800 


4, 


.00 


1, 


.19 








0. 


,022 


3900 


3900 


8. 


15 


0, 


,58 


0, 


.49 


0, 


.36 


VZ Lib 





.358 


0. 


,237 





.61 


5800 


5800 


3, 


.94 


1, 


.21 


0, 


.2 




0. 


,045 


4700 


4200 


7. 


,31 


0, 


,67 


0, 


.55 


0, 


.45 


V2388 Oph 





.802 


0. 


186 


0, 


.41 


6100 


6800 


1, 


.78 


2, 


.10 


0, 


.19t 


0.088 


0. 


10 


5900 


5900 


3. 


,59 


1, 


,36 


0, 


.65 


0, 


.55 


BB Peg 


0, 


.362 


0. 


,360 


0, 


.52 


5900 


6200 


3, 


.65 


1, 


.33 








0. 


,009 


4000 


3900 


8. 


,76 


0, 


,51 


0, 


.39 


0, 


.28 


HX UMa 





.379 


0. 


.291 


0, 


.44 


6600 


6700 


3. 


.32 


1, 


.44 


0, 


.047t 


0.63 


0. 


,023 


4400 


4400 


6. 


,64 


0. 


.73 





.50 


0. 


.39 


II UMa 





.825 


0. 


.172 


0, 


.4 


6600 


6800 


1. 


.70 


2, 


.15 


0. 


.23t 


0.87 


0. 


,15 


6100 


6400 


3. 


.29 


1. 


.45 





.67 


0. 


.58 


HT Vir 





.408 


0. 


.812 





.56 


6100 


6000 


3. 


.54 


1, 


.23 


0. 


.586§ 


0.6 


0. 


,28 


6100 


5900 


4. 


.12 


1. 


.15 





.93 


0. 


.52 


SV Camd 





.593 


0. 


.641 


0, 


.62 


5800 


5800 


3. 


.00 


1, 


.49 








0. 


.016 


3900 


3900 


7. 


.49 


0. 


.65 





.44 


0. 


.27 


SW Lynd 


0, 


.644 


0. 


.524 


0, 


.38 


7200 


6900 


2. 


.12 


1, 


.87 


0. 


.33 




0. 


.19 


6200 


6500 


3. 


.90 


1. 


.23 


0, 


.66 


0. 


.43 


KP Pes'* 





.727 


0. 


.322 





.06 


7400 


8900 


0, 


.92 


2, 


.68 


0. 


.52t 


3.5 


0. 


OS"^ 


7700 


6600 


1. 


.63 


2. 


.19 





.82 


0. 


.62 



''For resolved triples, was used to infer Tq and Mv\ for all others, /3fit was used. 

''The values inferred from our fit are inaccurate, since the tertiary is brighter than the contact binary and dominates the average spectrum. 
'^The flux ratio is inaccurate, since only a fraction of the tertiary's light fell inside the slit. 

''SV Cam, SW Lyn, and KP Peg are not contact binaries. Therefore, the components do not have equal temperature and the deduced properties are less 
reliable. They are not used in our statistical analysis. 



References. — The binary periods (P), mass ratios (M2/M1) and B — V are taken from the original papers I- IX (see Table 1). Put. and the separation 
angle are from the same sources except where indicated: ]: Tokovinin 1997; f: Paper I; §: (Heintz 1986). 



In order to verify the above, we also derived ter- 
tiary temperatures from My (column T/j). These 
should be similar to those inferred from the spec- 
tral fits (Tfit); from Table 2, one sees that this is 
indeed the case. We note, however, that this is 
not a strong test, since temperature does not vary 
strongly with stellar mass.^ 

Finally, in order to derive mass ratios, we also 
need an estimate for the masses of the stars in 
the contact binary. For this purpose, wc use that 
the contact binary's luminosity will be the sum of 
the luminosities of its two roughly main-sequence 
components. Since the luminosity depends steeply 
on mass, this implies that for low mass ratios 
q = M2/M1, the luminosity is simply that of the 
primary, while for higher ones there is a contri- 
bution from the secondary. To estimate the pri- 
mary's absolute magnitude, My.i, we use that for 
main-sequence stars, the V-band luminosity scales 
as M^-^ (inferred from Cox 2000, Table 15.7 and 
15.8; for a more detailed analysis, see Mochnacki 
1981), so that 

My,i ~My-K2.51og(l + f/-4) (2) 

Next, we estimate the primary's mass Mi using 
the main-sequence mass-luminosity relation (Cox 
2000, Tables 15.7 and 15.8). Including the uncer- 
tainty in how far the star has evolved on the main 
sequence, we expect these masses to be accurate 
to <20%. 

4. STATISTICAL ANALYSIS 

In order to determine whether the number of 

tcrtiaries we find is consistent with the hypothesis 
that all close binaries form in triple systems, we 
need to ensure our sample is homogeneous. Thus, 
we need to remove systems for which our proce- 
dure did not work properly, and consider for which 
separations and masses (or mass ratios) we can be 
certain we would have detected a tertiary if one 
were present. We flag all systems we exclude from 
our sample in Table 1 and discuss our reasons in 
more detail below. Next, wc compare our results 
for tertiaries with those found for secondaries for 
solar-type stars, trying to extrapolate the tertiary 
frequency to masses and separations to which we 

^For the same reason, it is not very useful to infer masses 
from the fitted temperatures. 



are not sensitive, and testing the hypothesis that 
all contact binaries are in triple systems. 

4.1. Limits and Biases in our Sample 

Among the sample of close binaries observed 
at DDO most are contact binaries, but nine are 
not: CN And (somewhat uncertain; Rucinski et al. 
2000), SV Cam, RZ Dra, SV Equ, FS Leo, SW 
Lyn, KP Peg, DV Psc and VI 130 Tau. For these, 
the assumption of a single spectral type for both 
stars is inappropriate and hence our procedure will 
not work optimally."^ We thus exclude all nine 
systems in our statistical analysis. 

We exclude a further five systems because the 
fit to the binary's spectrum was too poor to detect 
a third star down to flux ratios of 0.008. For one 
of these, V523 Cas, the temperature is low and we 
do not have a good template in our library. For 
a further four, AH Aur, DK Cyg, BX Dra, V351 
Peg, the match to the line profile is poor. 

We now turn to physical limits and biases. 
First, since our method is based on spectra taken 
through a 1'.'8 slit, wc will only be able to detect 
tertiaries at relatively close separations. For sepa- 
rations in excess of ~ I'.'S, the contribution of ter- 
tiary light will be reduced and hence wc will only 
detect very bright objects (such as V776 Cas). At 
a typical distance of ^ 100 pc, and including a sta- 
tistical correction factor of 10°'^^ = 1.35 (as in 
Duquennoy & Mayor 1991) for projection efi'ects, 
this corresponds to a separation of ~ 240AU or, 
assuming a total mass of the system of ^ 2Mq, 
an orbital period of ^2600yr ~ 10^ d. 

Second, a tertiary needs to be sufficiently 
bright. For main-sequence stars, we cannot de- 
tect tertiaries with flux ratios below ^ 0.008. 
From Table 2, one sees that this corresponds to 
mass ratios M3/M12 — 0.28. The value does 
not appear to depend much on the properties 
of the contact binary. To see why, we consider 
three possible configurations spanning the ex- 
tremes of the range of contact-binary proper- 
ties seen, with masses Mi = {1,1,2}M© and 
M2/M1 = {0.3,0.8,0.3}. For those parameters, 
the absolute magnitudes Mva2 — {4.7,4.4,2.0} 
and a tertiary with, e.g. /3—l% would be at 

^We beheve the detections of tertiaries in SV Cam, SW Lyn, 
and KP Peg are reliable despite the fact that these systems 
are not in contact. 
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Mv,z =i {9.7,9.4,7.0}. This corresponds to ter- 
tiary masses M3 ~ {0.42, 0.45, 0.70} M©, and thus 
to mass ratios M3/M12 ~ {0.32,0.25,0.27}. We 
conclude that we could detect all tertiaries with 
mass ratios in excess of 0.28. 

Third, independent of our method, we must 
somehow be able to observe a contact binary. For 
very bright tertiaries, the contact binary would be 
completely outshone and it likely will not be de- 
tected unless the tertiary is a star that appears 
interesting on its own accord and is studied in de- 
tail. But there is a bias even if the tertiary is 
less bright, contributing, say, only half the flux. 
In such a case, the variability would still be de- 
tectable, but the system might well be misclassi- 
fied, since the narrow lines of the tertiary would 
stand out in the spectrum while the broad ones 
from the contact binary would be much harder to 
detect (a good example of such a system is TU 
UMi; Rucinski et al. 2005). Given the above, we 
expect the sample of known contact binaries to be 
biased against systems having tertiaries brighter 
than the contact binary (which, roughly, corre- 
sponds to a mass in excess of that of the primary, 
or a mass ratio M3/M12 > 0.75). 

Finally, more generally, since we rely on the 
spectroscopic signature of a tertiary, we cannot 
detect white dwarfs, neutron stars or black holes. 

In summary, out of a sample of 75 close binaries, 
there are 61 contact binaries for which our method 
worked well, among which we detect 20 tertiaries. 
Among these, however, three (44 Boo, V899 Her 
and KR Com) should not be included, since the 
tertiary is at least half as bright as the binary and 
wc cannot be confident our sample of contact bi- 
naries is unbiased for such systems. Furthermore, 
V776 Cas and GZ And should not be counted as 
tertiaries, since the separation between the binary 
and third star is too large, making us incomplete. 
Thus, we are left with a sample of 58 contact bi- 
naries, for which we can be reasonably confident 
that our 15 detected tertiaries constitute all main- 
sequence tertiaries with 0.28 < Ms/Mcb $0.75 
and P3 < 10'^ d. 

4.2. Comparison with Solar-type Binaries 

Our method will miss triple systems with ter- 
tiaries that are at large separations, have low mass, 
and/or are compact. To estimate their number we 



need to extrapolate, but we do not know a priori 
the mass-ratio and separation distributions of the 
tertiary. By way of an estimate, we will assume 
that these are the same as those found for binary 
companions to solar- type stars. Along the way, we 
will try to test the alternate hypothesis that not 
all contact binaries have tertiaries, but that rather 
that the companion frequency is similar to that of 
solar-type stars. 

Before proceeding, we note that our choice of 
alternate hypothesis is somewhat arbitrary. One 
would like to test the hypothesis that the number 
of tertiaries we find is consistent with what one ex- 
pects from the formation of multiples. This num- 
ber, however, is not known: multiplicity among 
very young stars is very poorly constrained, and 
even among older stars there is no complete cen- 
sus (for a status report, see Tokovinin 2004). Our 
choice of comparing with solar-type stars corre- 
sponds to an implicit assumption that the com- 
panion frequency is independent of whether or not 
an inner system is a single star or a binary. It 
seems likely that this is a conservative assumption, 
i.e., assuming multiplicity plays no role in the for- 
mation of close binaries, the companion frequency 
of contact binaries is unlikely to be higher than 
that for solar-type stars. 

For the solar-type stars, we use the Duquen- 
noy & Mayor (1991, hereafter DM91) sample of 
164 stars of spectral type F7 to G9 (masses ^ 0.8 
to 1.3 M0). Using their distributions of mass ra- 
tio and period, we will mimic the selection effects 
present in our sample. Before doing so, however, 
two complications need to be mentioned. One is 
that among the 81 orbits used in the mass ra- 
tio and period distributions, six are second orbits 
from triple systems,^ and four are second and third 
orbits from two quadruple systems. For our pur- 
pose of estimating probabilities of finding a com- 
panion with certain parameters, this leads to an 
overestimate (e.g., for the full DM91 sample, the 
number of single stars is 93, not 164—81 — 83). To 
avoid biasing ourselves, we treat the DM91 sam- 
ple as consisting of 174 targets (some of which are 
stars, some binaries, some triples), 81 of which 
have a companion. 

A second complication is that the sample of 

^DM91 mention seven triple systems, but do not use the 
extremely wide outer orbit of the triple HD 122660. 
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DM91 is divided into two groups: spectroscopic 
and resolved (visual and common proper motion) 
binaries, which correspond to binaries with pe- 
riods P < 10^ d and P > 10"* d, respectively. 
These groups differ in that white dwarfs will be 
present among the (single-lined) spectroscopic bi- 
naries, but not among the resolved ones. In their 
statistical analysis of mass ratios, DM91 inserted 
eight 'fake' white-dwarf systems in the long-period 
group: two each in the four mass bins with (/max = 
{0.5,0.6,0.7,0.8}. We will remove these, since we 
cannot detect white dwarf companions. We will 
also correct for the presence of white dwarfs among 
the spectroscopic binaries. 

We now turn to the application of the selection 
effects present in our sample. First, to mimic the 
incompleteness among known contact binaries, we 
ignore all binaries with mass ratios q = M2/M1 > 
0.75: from lines (3) and (5) in Table 7 of DM91, 
after correction for one 'fake' white dwarf in line 
(5), we find this reduces the sample by 6.15 spec- 
troscopic and 7.5 resolved binaries. (Here, the 
numbers are non-integer since we had to split the 
q = 0.7-0.8 bin and since for the spectroscopic bi- 
naries, DM91 corrected for the distribution of or- 
bital inclinations). This leaves a sample of 160.35 
targets. 

Second, to reproduce our detection limit, we 
count all binaries with mass ratio q > 0.28 among 
the remaining targets: again from lines (3) and 
(5), after correction for seven 'fake' white dwarfs, 
we find 22.27 and 31.8 binaries, respectively. 

Third, to account for our separation limit, we 
select systems with periods P < 10^ d. For this 
purpose, we use that from the period distribu- 
tion (Fig. 7 in DM91), among the resolved sys- 
tems with P > 10^ d, 31 out of 65 have P < 10^ d. 
Thus, only 15.17 out of the 31.8 long-period bi- 
naries remain, and the total implied detection 
rate for a survey like ours would be (22.27 -|- 
15.17)/160.35= 23%. 

In the above, we still need to correct for the 
presence of white dwarfs among the spectroscopic 
binaries. DM91 mention that from statistical 
considerations of stellar populations, one expects 
"about two white dwarfs per decade of period." 
This would imply about eight are present among 
their sample of 34 spectroscopic binaries. Seven of 
these would be included in the 22.27 binaries with 
0.28 < q < 0.75 selected above, implying a re- 



duced detection rate of 19%. This may be an over- 
estimate. On the other hand, the selected DM91 
sample includes a number of companions with pe- 
riods below 10 d, which could not exist around con- 
tact binaries.^ Since these numbers are no more 
than guesses, we will use a rounded expected de- 
tection rate of 20% for our analysis below. 

In summary, we conclude that if companions to 
contact binaries occurred at the same frequency 
as those to solar-type stars, and if their properties 
followed the same mass-ratio and period distribu- 
tions, we should have detected tertiaries for about 
20% of our sample, or 12 out of 58. In reality, we 
found 15 tertiaries. This is 3 more than expected, 
but the difference is not highly significant: there 
is a 17% probability to find 15 or more tertiaries 
out of 60 systems when the expected tertiary rate 
is 20%. 

To calculate a similar probability for our hy- 
pothesis that all close binaries are in triple sys- 
tems, we need to estimate the probability that 
a tertiary will have the correct properties to be 
detected with our method. For this purpose, we 
use again the DM91 sample and first estimate the 
companion fraction for all systems with q < 0.75, 
again by adding up lines (3) and (5) in Table 7 
(we thus include the white dwarfs). Without the 
highly uncertain q < 0.1 bin, we find 34.75 spec- 
troscopic and 59 resolved binaries, respectively, or 
an implied companion fraction of 93.75/160.35 = 
58%. If we include the estimated 5.6 spectroscopic 
and 14 resolved binaries with q < 0.1, this rises 
to 67%. Thus, between 30 and 34% of all com- 
panions have q > 0.28, P < 10^ d, and are not 
white dwarfs. If the tertiary rate were 100%, these 
would be the expected detection rates, and hence 
we would expect to have found between 17 and 20 
tertiaries in our sample of 58. Conservatively as- 
suming the expected tertiary rate is 34%, we find 
that there is a 12% probability of finding 15 or 
fewer systems. Thus, our results are also consis- 
tent with the hypothesis that all close binaries are 
in triple systems. 

Finally, for all our above estimates we as- 
sumed that the tertiaries followed the same 
mass-ratio and period distributions as those of 
solar-type companions. We do not have suf- 
ficient objects to test this rigorously, but can 

® One of these, 44 Boo, is in our sample as well. 
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at least verify this hypothesis. For the ranges 
q = {[0.28, 0.4), [0.4, 0.5), [0.5, 0.6), [0.6, 0.75)}, we 
found {8,3,3,1} tertiaries. Consulting Table 7 
of DM91, and scaling to the same total num- 
ber of companions (15), we infer {6.4,3.2,3.5, 1.9} 
binaries (where, as above, we reduced the long- 
period bin by a factor 31/65 to correct for peri- 
ods > 10^ d and deducted one system in the three 
higher mass-ratio bins in order to correct for white 
dwarfs among the spectroscopic binaries) . Clearly, 
within the limited statistics, the two distributions 
are consistent. 

5. CONCLUSIONS AND PROSPECTS 

We searched a sample of 75 close binaries for the 
spectroscopic signature of tertiaries and identified 
23 triple systems, implying a ratio of almost one 
in three. For a homogeneous subset of 58 contact 
binaries, we are fairly confident our 15 tertiaries 
are all those that have periods < 10^ d and mass 
ratios 0.28 < Mg/McB < 0.75. 

We compared our results with expectations un- 
der two hypotheses, that the incidence of tertiaries 
is similar to the incidence of companions to solar- 
type stars, and that all close binaries are in triple 
systems. The latter hypothesis is expected to hold 
if close binaries form via the Kozai mechanism; 
the former appears a conservative upper limit for 
the case that the formation of close binaries is un- 
related to multiplicity. Using the DM91 sample 
of companions to solar-type stars to infer mass- 
ratio and period distributions, we find that, for 
the two hypotheses, the expected numbers of triple 
systems among our sample are 12 and 20, respec- 
tively. Finding 15 systems is consistent with either 
hypothesis. 

While inconclusive in terms of testing the role of 
multiplicity in the formation of close binaries, we 
feel the relatively large fraction of triple systems 
found is encouraging, especially since in Paper 1, 
from a variety of methods, a high tertiary fraction 
of 59 ±8% was observed as well. To make progress, 
a larger fraction of tertiary parameter space needs 
to be covered. For the method presented here, this 
is not difficult since the archive observations used 
here were not optimized for the search for faint ter- 
tiaries. By using higher resolution spectra span- 
ning a larger wavelength range, one can improve 
the contrast in the spectra and increase the sensi- 



tivity, and by observing at longer wavelength, one 
will be sensitive to lower mass tertiaries for a given 
limiting contrast ratio. All three improvements 
are possible with echelle spectrographs. Further- 
more, with adaptive optics in the near infrared, 
one can reach even lower mass tertiaries (though 
only on relatively long orbits). We hope to follow 
both routes in the future. 

We thank Peter Eggleton and Andrei Tokovinin 
for enlightening discussions, and Stefan Moehnacki 
and Fang Bao for help in the initial stages of this 
project. We made extensive use of the SIMBAD 
database and the VizieR catalogue access tool, 
both operated at CDS, Strasbourg, France, and of 
NASA's Astrophysics Data System. We acknowl- 
edge financial support by NSERC. 
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A. INDIVIDUAL SYSTEMS 



Below, we briefly summarize; tiic; properties of all binaries for which wc; detected the spectroscopic signature 
of a tertiary, as well as for a number of systems without detections for which the results warrant further 
discussion (marked with f in flux ratio column, /3, of Table 1). The roman numeral directly following the 
system's name is the paper in the series of "Radial- Velocity Studies of Close Binary Stars" from which our 
data were drawn; see the references at the end of Table 1. For data on resolved multiples, we generally rely 
on the Multiple Star Catalogue (MSC, June 2005 update; Tokovinin 1997). 

GZ And (I) is a W UMa binary, and the brightest component of a visual multiple system. Despite the 
noisiness of the spectrum and the relatively poor match for our best temperature of 5600K, we clearly detect 
a third companion, with /3 = 0.015 and T3 = 3500 K. Of all visual components, the only one that this could 
correspond to, is component E of Paper I, which is at a separation of 2'.'13 and has AH = 2.6, AK = 2.4. 
Using Table 2 and Cox (2000, Tables 7.6, 15.7, and 15.8), we infer Mk = 3.3 for GZ And. Thus, component 
E has Mk ^ 5.7 and is likely an ~O.42M0 M2 star with T3 ~ 3500 K, V - K ~ A.O and Mv,3 9.7. The 
temperature agrees very well with our measurement, and hence we are confident we detected component E. 
The implied flux ratio of /? ~ 0.011 is somewhat smaller than what we measure, the opposite of what is 
expected given that the tertiary would not have been completely in the slit. Likely, our measurement is 
biased by the relatively poor fit. 

HV Aqr (III) is an A-type contact binary system, and is one of the best examples of our program's ability 
to detect components at fiux ratios of only a few percent. Although visually there appears little improvement 
between the fit with and without a third star and it is difficult to see the third star's contribution in the 
residuals from the binary-only fit, there are distinct minima in the variance as a function of T3 and Av^. 
Furthermore, the fitted temperature, T3 [n = 4000 K, matches that inferred from the flux ratio, T^^p = 4000 K. 
In Paper I, no other indicators for multiplicity were found. 

V410 Aur (VIII) is a W UMa-type binary in a known triple system, with a tertiary at 1'.'7 that is fainter 
by AV = 1.04 (MSC), corresponding to /3 = 0.38. The signal of the tertiary is obvious in the spectra and 
from its contribution to the broadening function, a flux ratio of 0.26 was inferred, while from our routine we 
infer 0.22. Likely, both numbers are lower than the true flux ratio because some of the light fell outside the 
slit. The temperature T^^m = 5200 K is consistent with that inferred from the flux ratio, T^^^ = 5600 K. 

44 Boo B (IV) is the contact binary nearest to Earth, and its spectrum is dominated by a brighter star 
{AV = 0.78) at a separation of 1'.'7. In the spectra used to analyze this system, some of the light from the 
third component was blocked by the slit, resulting in a flux ratio of 0.4-0.7 inferred from the broadening 
function. Our procedure yields a somewhat lower flux ratio of 0.23, although the program fits part of the 
third star as if it were the contact binary. The fit yields very sharp, well-defined minima in variance as a 
function of T3 and Av^, but our results are nevertheless poorly defined, since the initial fit to the contact 
binary was biased greatly by the presence of the third star (since it dominates the spectrum). As a result, 
the temperature inferred for the contact binary is too high, and that for the tertiary too low. Since we use 
the observed fiux ratio, however, our inferred tertiary mass and mass ratio should be accurate. 

CK Boo (II) is an A-type W Uma system. We find a good fit to the spectrum, and while there are some 
systematic residuals, a clear signatiire of a third star is present, with a very low flux ratio /? = 0.009. The 
fitted temperature is consistent with that inferred from the flux ratio, T^^/j = 3900 K. The tertiary is also 
detected in the adaptive optics observations described in Paper I. The separation is only 0'.'12 and hence 
the magnitude difference AK ~ 2.8 is rather uncertain. We nevertheless tried to verify consistency: using 
Table 2 and Cox (2000, Tables 7.6, 15.7, and 15.8), we infer Mk = 2.6 for CK Boo and thus Mk,3 ^ 5.4. 
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The latter implies that the tertiary would be an -O.48Af0 Ml star with T3 ~ 3800 K, V - /-C ~ 3.8 and 
Mv,3 — 9.2. The temperature and implied flux ratio of 0.007 agree well with our measurements. We note 
that from arrival-time variations, in Paper I the possibility of a companion in an ^5 AU orbit was mentioned. 
The inferred minimum mass of 1.5 M©, however, is inconsistent with our results, unless it were a neutron 
star or black hole. 

FI Boo (IV) is a W-type contact binary system. We clearly detect a faint third component, with /3 = 0.012 
and T3 (it = 3900 K. The fit to the binary has relatively large systematic residuals, which dominate the 
variance; as a result adding a third star does not change the variance as much as might be expected to if the 
fit were better. Despite these limitations, there is an obvious minimum in of the variance as a function of T3 
at a much later spectral type than that of the main binary. Furthermore, the variance shows a sharp drop 
at Av3 ~ 0. The presence of a tertiary is also inferred from stochastic residuals in Hipparcos measurements 
(Paper I). 

S V Cam (VI) is a detached binary, for which the likely presence of a third body in a 41 or 58 yr orbit was 
inferred from arrival-time variations (Lehmann et al. 2002; Borkovits et al. 2004); the implied separation is 
a few O'.'l and the mininiimi mass is around 0.2 Mq. The different temperatures of the binary components 
make it less suited to our method of analysis, and our fit to the average spectrum is relatively poor. We 
nevertheless clearly detect a tertiary, and infer a temperature Ta^^ = 3900 K that agrees nicely with the 
one derived from the fit, T3 gt = 3900 K. We note, however, that because of the poor fit, there is a marked 
decrease in variance for all tertiary spectral types: likely, the third star is being fit to some of the residuals 
left from fitting the main binary. As a result, the fitted flux ratio of 0.016 may be somewhat higher than it 
would be if it were fitting only the third star. Since the temperature is close to the lower limit of our range 
of templates, the real temperature may well be lower. An independent indication for a lower temperature 
would be that the inferred mass of M3 = 0.65 Mq is somewhat high compared with that inferred from the 
arrival-time orbit: it would require an inclination ^3 < 20°, which has an a priori probability of < 5%. 
Furthermore, Lehmann et al. (2002) mention that "masses of >O.6OM0 should be excluded because a third 
stellar spectrum would be visible in the observations which is not the case." Since SV Cam is not a contact 
binary, we do not include it in our statistical analysis. 

V776 Cas (V), an A- type contact binary, is the brighter member of a visual binary, with an angular 
separation of 5'.'38 and AV = 1.56 (MSG). Our fit to the contact binary is fair, and easily detect the third 
star. Our flux ratio is much lower than the measured one since most of the tcrtiary's light fell outside of 
the slit (indeed, many of the individual spectra of V776 Cas were taken on purpose excluding the third star 
to make it easier to calculate the radial velocity for the individual components). The fitted temperature 
Tafit = 6100 K is somewhat higher than that inferred from the flux ratio. 

EE Get (VI) is a component of a visual binary in which the second component, at 5'.'6, is another close 
binary (hence, the system is a quadruple). The contact binary is the fainter component, by AV^ = 0.36 
(MSG). The data used in our observations were taken excluding as much light from the other binary as 
possible, but our program nevertheless picks out scattered light from the companions and identifies it as a 
very faint third star. We do not list it in Table 2, since without knowing that the system was a multiple, wc 
would not have identified it as such. For reference, we note that our fit yields T12 = 6200 K and T3 = 6800 K. 

V445 Cep (IX) is an A-type contact binary with very shallow eclipses. It is a hot system - we measure 
T12 = 7400 K for the main binary. Wc also find a third star in the system, with a T3 = 5600 K and relative 
flux (3 = 0.055. This flux is high enough that one would have expected the tertiary to have been detected in 
previous surveys, so it is a bit puzzling that it has not. Nonetheless, the fact that there is a clear minimum 
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in the variance wlien a third star of a very different spectral type is added leads us to conclude that there is 
a third component in the system. 

KR Com (VI) is an A-type contact binary in a visual binary, with a companion at C/.'lld that is fainter 

by AV — 0.59, or /? = 0.58 (MSG). From the broadening function, a flux ratio of 0.56 was inferred, while 
we find a much lower value of 0.23, likely because our fit is biased by the fact that the tertiary is so bright. 
Because of this, we use the observed flux ratio to infer the tertiary's parameters. 

V401 Cyg (VI) is a contact binary for which the spectral signature of the tertiary, despite being only 
at the 3% level, was already seen in the broadening function. Our fit to the binary is poor in a somewhat 
surprising fashion: in some parts, it reproduces the spectrum very well, while in others, particularly around 
5170 A, it fails utterly. Despite the resulting uncertainties, we clearly recover the third star with /3 = 0.015. 
The temperatures inferred from the fit and the flux ratio are both 4700K. The presence of a close-in tertiary 
is also inferred from stochastic residuals in Hipparcos measurements (Paper I); it cannot be the object at a 
separation of 18'.'0 in the adaptive optics observations presented in Paper I. 

V2082 Cyg (IX) is likely an A-type contact binary, although a detached configuration cannot be com- 
pletely excluded. We measure T12 = 7000 K, and obtain a fairly good fit, although there are a number of 
features - at 5235, 5195, and 5167A - that are stronger than in the template. Nevertheless, we are able to 
detect a faint third companion, with /3 = 0.02 and T^^t = 5100 K (the latter consistent with T^^fj = 5200 K 
inferred from the flux ratio) . The only other indication for the presence of a tertiary found in Paper I, was 
that the X-ray flux was stronger than expected for the early-type binary. 

V2150 Cyg (IV) is an A-type contact binary that has a much fainter, AV = 3.35 visual companion 

at 3'.'68 (Paper I). The faint companion star was outside the slit in our spectra, and unlike for brighter, 
well-separated visual doubles in our sample (such as EE Get), we detect no scattered light from the third 
companion, nor do we find any closer companion. Our detection limit is at the 1% level, despite the fact that 
the binary's spectrum is not fit all that well, with relatively large, broad-scale deviations as well as some 
higher-frequency residuals. This may partly be because the temperature of about 7200 K places V2150 Gyg 
near the upper end of our temperature range, where we have relatively few template spectra. 

V899 Her B (IV) is an A-type contact binary for which the broadening profiles indicate it is the fainter 
component of a spectroscopic triple: (3 = 1.5. The brighter star is itself a radial-velocity variable as well. 
Prom our procedure, we infer a smaller value for the flux ratio. 0.725, but this is likely because the initial fit 
to the spectrum is biased by the dominating flux from the third star. Hence, the temperature for the binary 
is overestimated. Nevertheless, the overall fit is fairly good. In Paper I, no other evidence for the presence 
of a tertiary is listed. 

ET Leo (VI) is a low-amplitude contact binary presumably seen at low inclination. Our initial fit to the 
binary spectrum left rather strong, large-scale residuals; these were reduced but not altogether removed 
using a tenth-order polynomial fit to the continuum. Our best-fit temperature for the binary is 5800 K. This 
is substantially higher than what would be inferred from the spectral type of G8 assigned by Rucinski et al. 
(2002), but consistent with the temperature inferred from B — V. The residuals from the best fit show a 
clear signature of an M-type tertiary, and this is confirmed by the variance as a function of T3 and AU3. The 
inferred fiux ratio and temperature are (3 = 0.022 and T^^i = 3900 K. The presence of the tertiary was also 
suspected based on residuals of Hipparcos measurements (Paper I). 
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VZ Lib (IV) is a contact binary for which the presence of a fainter tertiary component was already 
indicated by its signature in the broadening function. No other indicators for a tertiary were found in Paper 
I. We recover the third star unambiguously, but find /? = 0.045, which is much fainter than the value of 0.2 
estimated from the broadening function. This may be because the previous measurement effectively measured 
the tertiary's flux assuming it had the same spectral type as the binary, while in reality it is later and hence 
has stronger lines (Sect. 3.1). Our overall fit is good, although some of the sharp-line features (particularly 
at 5182 A) are reproduced relatively poorly, indicating that the tertiary spectral type or metallicity may not 
be entirely correct. The tertiary's temperature estimates also not quite in agreement, T^^/} = 4200 K from 
the flux ratio and Taju = 4700 K from the temperature. 

SW Lyn (IV) is a close, semi-detached binary in a 'reversed-Algol' configuration, with the more massive 
component filling its Roche lobe. Despite the fact that the components are not in contact and therefore have 
different temperatures, we easily recover the tertiary that was identified from the broadening function. We 
find /? = 0.19, substantially less than the value of 0.33 inferred before (likely because the latter does not take 
into account that the tertiary has a different spectral type; Sect. 3.1). Overall, the fit is good, although some 
mismatches remain. Since this system is not a contact binary, it was not used in the statistical analysis. 

V502 Oph (IX) is a W-type contact binary, for which there are several pieces of evidence pointing to a 
third companion. First, Hughes & McLean (1984) detected two radio sources near the source, separated 
by only 2'.'6. Second, Derman & Demircan (1992) found that the arrival times of the minima showed a 
modulation with a period of about 35 years (which, however, is too short for a companion separated by 
2f.'6). Third, Hendry & Mochnacki (1998) detected stationary Na I lines in trailed spectra. The latter 
detection is the most convincing, and is the basis of the identification as a triple in Paper I. Unexpectedly, 
our procedure does not unambiguously recover the tertiary. We find a reasonable fit for a binary temperature 
Ti2,fit = 5800 K, but with some rather odd residuals throughout the spectrum, indicating our template is not 
a good match. Adding a third star, we do see an improvement in the quality of the fit, but the minimum in 
variance does not occur close to zero relative velocity. Furthermore, taken at face value, the relative fiux of 
/? = 0.007 indicates a tertiary temperature much lower than the fitted value, Ta = 6900 K. For this reason, 
we have not counted this system as a detection. We note, however, that the discrepancy might be reduced if 
the separation is really 2'.'6, since in that case much of the tertiary's light might have fallen outside the slit. 

V2388 Oph (VI) is a W UMa member of a bright visual binary, with a separation of 0'.'087 and magnitude 

difference AV = 1.75 (MSG). It is possibly seen in arrival-time and astrometric variations as well (Paper I). 
From the broadening function, a tertiary fiux ratio /3 = 0.2 was found, while our procedure yields a value 
of 0.10. We find we can reproduce the spectrum very well, although the initial fit (without a third star) 
largely incorporates the light from the third star. Hence, the inferred binary temperature, Ti2jii = 6100 K, 
will be biased somewhat (the fact that the temperature inferred from B — V is similar likely refiects the fact 
that the colour is contaminated by the tertiary as well). Nonetheless, we find a very different spectral type 
for the tertiary, with T^^^ = 5900 K, consistent with what is inferred from the flux ratio. 

BB Peg (I) is a W-type contact binary in which the combination of a light-curve fit and radial-velocity 
orbits allowed the masses of both components to be measured: Mi = 1.38 and M2 = 0.5 Mq. The 
primary mass is in good agreement with the mass inferred from the absolute magnitude through the period- 
himinosity-color relation. The average spectrum, which is among the more noisy we analyzed, is reproduced 
fairly well by our procedure, although some small deviations on relatively large scales remain. There is a 
clear drop in the variance as cooler third stars are added to the system, indicating the presence of a third 
component in the system with a relative fiux /3 = 0.009 and temperature T^^^ = 3900 K. A third component 
is also suspected from arrival-time variations (Paper I). 



19 



KP Peg (IX) is a /3-Lyrac type binary, and is the brighter component of a visual binary with a separation 
of 3'.'5 and AV = 1.6. It has one of the earliest spectral types, A2, implying a temperature that is outside 
the range covered by our templates. Nevertheless, we find a fairly good fit to the binary, and easily detect 
the third star. We find [3 = 0.031, which is dimmer than the known flux ratio since most of the light of the 
third star did not enter the slit. The fitted tertiary temperature of T^^fit = 7700 K is not consistent with 
what is expected for the observed magnitude difference. 

V335 Peg (IX) is an A-type contact binary with the second component contributing only 5% of the total 
flux. We measure T^\2 — 6400 K, and see some minor residuals. The variance profile for the third star 
declines for later type stars, indicating a third star with T3 = 3700 K. The estimated flux it finds for the 
third star is only (i =0.006, however, which is right at our detection limit. Hence, we classify it as an 
interesting non-detection. 

HX UMa (VIII) is a A-type contact binary with a previously identified, fainter (AF = 3.31) companion 
at a separation of 0'.'626 (MSG). The tertiary 's signal was also detected in the broadening function, yielding 
(3 = 0.049. Our procedure provides an excellent fit to the main binary, with the fit improving even further 
with the addition of a third star. We find (3 = 0.023, with is somewhat lower than indicated by the magnitude 
difference and inconsistent with the measurement from the broadening function. It is possible that this results 
partly from light not entering the slit and partly from the bandpass being blueward of V. (If so, the good 
agreement foimd earlier would be due to a fortuitous cancellation of the light loss by the overestimate of the 
flux resulting from the assumption that the tertiary had the same spectral type as the binary; Sect. 3.1.) 
The fltted tertiary temperature is T^^^ = 4400 K, which is in excellent agreement with that inferred from 
the observed flux ratio. 

II UMa (VI) is an A-type contact binary that has a fainter, AV = 1.64, companion at a separation of 0'.'87 
(MSG). The tertiary was obvious already in the broadening function, yielding /3 = 0.17, and our procedure 

easily recovers it. We find f3 = 0.15. Unlike other cases, this is consistent with the value found from the 
broadening function, since the spectral types of the contact binary and the tertiary are very similar. 

HT Vir B (IV) is part of a close visual binary, with a period of 274 yr and semi-major axis of I'.'Ol (Heintz 

1986), current separation ^ 0'.'6), and the spectrum shows strong lines from the third star. The contact 
binary is brighter than the companion during its maxima (AV^ = 0.63), but fainter during the minima. Lu 
et al. (2001) found radial-velocity variations in the third component, indicating that it is in a binary itself 
and hence that the system is a quadruple. Our procedure yields a flux ratio f3 = 0.28, lower than the value of 
0.52 inferred from the broadening function. Our fit is fairly good, and we infer temperatures around 6000 K 
for both components, consistent with the values inferred from the colours and flux ratio. 
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